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We determine the spectrum of scalar masses in a supersymmetric vacuum of a general 
^ J\f = 1 supergravity theory, with the Kahler potential and superpotential taken to be random 

^ functions of N complex scalar fields. We derive a random matrix model for the Hessian matrix 

and compute the eigenvalue spectrum. Tachyons consistent with the Breitenlohner-Preedman 
* bound are generically present, and although these tachyons cannot destabilize the supersym- 

metric vacuum, they do influence the likelihood of the existence of an 'uplift' to a metastable 
^ (— I vacuum with positive cosmological constant. We show that the probability that a supersym- 

metric AdS vacuum has no tachyons is formally equivalent to the probability of a large fluc- 
^ tuation of the smallest eigenvalue of a certain real Wishart matrix. For normally-distributed 

^ matrix entries and any N, this probability is given exactly by P = e^p{—2N^\W\'^ /mly^gy), 

^ with W denoting the superpotential and mgusy the supersymmetric mass scale; for more gen- 

eral distributions of the entries, our result is accurate when A S> 1. We conclude that for 
t — ■ 1^1 ^ f^susy/N , tachyonic instabilities are ubiquitous in configurations obtained by uplifting 

^ supersymmetric vacua. 
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1 Introduction 



Supersymmetric vacua of four-dimensional supergravity provide a computable stepping 
stone toward realistic non-supersymmetric vacua of string theory. Supersymmetry se- 
cures vacuum stability, and also opens powerful analytic perspectives, such as the 
AdSi/CFT^ correspondence, that are unavailable in de Sitter solutions. AdS vacua 
have played a prominent role in the development of mechanisms for moduli stabiliza- 
tion, as exemplified by [1]: a supersymmetric AdS vacuum serves as a controllable 
foundation, and a source of supersymmetry breaking is introduced to lift the energy 
density to positive values. However, tachyons are abundant in supersymmetric AdS 
vacua, and may pose significant obstacles to uplifting AdS solutions to metastable de 
Sitter vacua. 

In this work we use random matrix theory to obtain the scalar mass spectrum in a 
supersymmetric vacuum of a generic four-dimensional = 1 supergravity. We take the 
Kahler potential K and superpotential W to be random functions of N complex scalar 
fields, in a sense made precise in [2] . The Hessian matrix T-iis a. particular combination 
of the derivatives of W and and the techniques of random matrix theory can be 
used to determine the spectrum of eigenvalues of i.e. the scalar mass spectrum. 
Strong correlations among the eigenvalues make large fluctuations — for example, a 
fluctuation resulting in a tachyon-free spectrum — extremely unlikely at large N . 

The key results of this work are statements about the prevalence of tachyons in 
supersymmetric AdS vacua. AdS supersymmetry (see [3] for a review) ensures that 
the real and imaginary parts of each complex scalar field have distinct masses [4] , with 
splitting proportional to the AdS scale, set by \W\. The mass spectrum therefore 
consists of two branches. We show that the spectrum of each branch is determined 
by the eigenvalue spectrum of a particular ensemble of real Wishart matrices, and we 
give an analytic result for the full spectrum. The mass spectrum critically depends 
on the ratio of the AdS scale to the supersymmetric mass scale, denoted nisusy'- for 
^susy/\W\ — 7- oo all scalars have positive mass-squared, with spectrum given by the 
Marcenko-Pastur law [5], while for msusy/\W\ — >■ tachyons are endemic and the 
spectrum asymptotes to that of the Altland-Zirnbauer CJ ensemble [6]. 

Using a result of Edelman [7], we obtain the probability P of a fluctuation that 
renders all scalars non-tachyonic, in terms of the probability of a corresponding fluc- 
tuation in the associated Wishart ensemble: we find that P = exp{—2N'^\W\'^ / ml^^y) . 
Then, building on results from the Coulomb gas formulation of random matrix theory, 
we obtain the mass spectrum resulting from such a fluctuation. We conclude that AdS 
vacua with Tn?susy A^^|W^P are overwhelmingly likely to contain tachyons allowed 
by the Breitenlohner-Freedman (BF) bound. Finally, we determine the probability 



1 



that a de Sitter critical point obtained by uplifting a supersymmetric AdS vacuum is 
metastable, for a range of uplifting scenarios (see also [8]). 

The organization of this paper is as follows. In §2 we present the structure of the 
Hessian matrix in a supersymmetric AdS vacuum. After briefly reviewing the Wishart 
ensemble [9] and the Altland-Zirnbauer CI ensemble [6], we show that the eigenvalue 
spectrum of the Hessian matrix is determined by the spectrum of the Altland-Zirnbauer 
CI ensemble, which in turn can be related to the spectrum of a particular ensemble 
of Wishart matrices. In §3 we compute the spectrum of the Hessian matrix and de- 
termine the probability of fluctuations of its smallest eigenvalue. We then study the 
distribution of tachyon-free supersymmetric vacua and argue that a fraction oc of 
the parameter space contains no tachyons at all. In §4 we describe the implications 
of our results, focusing on instabilities in uplifted vacua. We conclude in §5. In the 
appendix we present the results of extensive numerical simulations that confirm our 
analytic findings. 

Unless otherwise specified, we work in natural units with Mp^ = StcGn = 1. 

2 The Hessian Matrix in a Supersymmetric Vacuum 

We begin in §2.1 by presenting the structure of the Hessian matrix in a supersymmetric 
vacuum. Then, in §2.2, we briefly recall two relevant ensembles of random matrices, the 
Wishart ensemble and the Altland-Zirnbauer CI ensemble, and we express the eigen- 
values of the Hessian matrix in terms of those of the Wishart and Altland-Zirnbauer 
CI ensembles. 

2.1 Supersymmetric vacua in A/^ = 1 supergravity 

Supersymmetric points satisfying Fa = DaW = {da + daK)W = are critical points 
of the J\f = 1 supergravity F-term potential, Vp = [FaF'^ — 3|H^p). The Hessian 
matrix at such a point is given by 

where 

/ Za' Z-,, -ZabW 

^ V -z-aiW Z-: Z,, 

and Zah = Zi,a = T^aDhW is the supersymmetric fermion mass matrix. Here Va denotes 
the geometrically and Kahler covariant derivative, and we have specialized to a Kahler 
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gauge in which {K) = at the critical point. The effect of the contribution proportional 
to the unit matrix is to shift the entire spectrum to more negative values by an amount 
2|VFp. Therefore, finding the spectrum of the Hessian in a supersymmetric vacuum 
amounts to determining the spectrum of the matrix l-Lzi which will be the focus of 
much of this work. 

To simplify T-Lzi we first write Z = U'EU'^, where f/ is a unitary matrix whose 
columns are orthonormal eigenvectors of ZZ, S = diag(Ai, . . . Xn), and the \a are real 
and nonnegative, with the eigenvalues of ZZ [10]. This Takagi factorization exists 
for any complex symmetric matrix Z. Performing the 2N x 2N unitary transformation 



Hz^U^HzU with U 



U 

u* 



(2.3) 



the matrix l-Lz can be written as 

'Hz 



(2.4) 



After rearranging the rows and columns in an obvious way, Hz takes the block diagonal 
form 



V-2 



( \\ -FAi 











-W\2 A^ 
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V 



(2.5) 



The spectrum of each block of 'Kz is 

a;,± = A^±|l^|A,. 



(2.6) 



Thus, the scalar mass spectrum in a supersymmetric vacuum is fully determined by 
the vev of the superpotential, W ^ and by the spectrum of eigenvalues A^ of ZZ. 

The eigenvalue uja- is minimized by \a = \W\/2, which gives Ua- = — |Vrp/4. So 
the smallest possible eigenvalue of Ti is 



\W\ 



-2|iyp 
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W\ 



3 

'i' 



-^\V\ = rngp, 



(2.7) 



where we have used that V = —3\W\'^ at a supersymmetric minimum, and mep denotes 
the Breitenlohner-Freedman bound [11, 12]. 
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2.2 Classical ensembles 

In §3, we will obtain the spectrum of the Hessian matrix using the methods of ran- 
dom matrix theory. In this section, we review the essential properties of the two 
matrix ensembles that are relevant for the analysis of §3: the Wishart ensemble and 
the Altland-Zirnbauer C/ ensemble. 

2.2.1 The Wishart ensemble 

The ensemble of Wishart matrices W [9] takes the form 

W = AA\ (2.8) 

where A is an x M real or complex matrix whose entries are independent and 
identically distributed (i.i.d.) random variables drawn from a statistical distribution 
with mean zero and variance a^, which we denote by fl{0, a). The probability density 
of the eigenvalues of W depends on the ratio M/N, and for our purposes it will suffice 
to consider the case M > N. 

Since a Wishart matrix is the Hermitian square of another matrix, it is necessarily 
positive semidefinite. Upon changing variables to an eigenbasis of W, the Jacobian 
of the transformation induces 'interaction terms' between the different eigenvalues ^a- 
For entries drawn from a normal distribution,^ the joint probability density is (cf. [18]) 



f P f 1 ^ AT N 

/(/ii,. . . ,/XAr) =C exp I -- I ^^/ia - 2 ln|/Xa - J^ln/ia 

\ \ a=l a<b a=l y 



, (2.9) 



where ^ = M — N + 1 — 2/P, and /3 = 1, 2 for real and complex matrices, respectively. 
In the Coulomb gas picture in which the joint probability density is interpreted as 
an exponential of the free energy of an ensemble of interacting particles, the non- 
negativity of a Wishart matrix corresponds to the presence of a hard wall at /i = 0. 

The probability density function (pdf) for the eigenvalues of W is given by the 
Marcenko-Pastur law [5], 

PMP (/^) = Na^j^ ^ -f^)if^-V-), (2-10) 

where 

r]±=Na\l±^)', (2.11) 

and 7] = M/N > 1. 



^Universality implies that more general distributions yield equivalent results at sufficiently large N 
[13-17], which we have verified directly in our simulations. 
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Figure 1. The eigenvalue spectrum pMp{f^) for an almost square Wishart matrix with 
M = Af + l = 6andfT = l/V^. 

In §3, we will be led to focus on a specific ensemble of Wishart matrices in which 
M = + 1, so that A is 'almost square'. By an abuse of language, we will likewise 
refer to the associated Wishart matrices W as 'almost square', though of course W is 
square for any M, A^. For this ensemble, rj = 1 + 1/N, so that to lowest order in a 
expansion, with a = 1/\/N, we have rj^ ^ and rj^ ^ A. A plot of the eigenvalue 
spectrum for this almost square case is shown in Figure 1. 

2.2.2 The Altland-Zirnbauer CI ensemble 

The Altland-Zirnbauer CI ensemble [6] consists of matrices of the form 



where Zab is a complex symmetric matrix with independent entries distributed as 



Equivalently, the Altland-Zirnbauer CI ensemble can be defined as the matrix ensemble 
satisfying 




(2.12) 



Zab G ^(0, 0") for a 7^ 6 and 



Zaa G ^(0, \/2a) (no sum on a) . 



(2.13) 



TMT^ = M* , 
CMC^ = -M , 
PMP^ = -M , 



(2.14) 
(2.15) 
(2.16) 
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Figure 2. The eigenvalue spectrum of the Altland-Zirnbauer CI ensemble for N = 100. 

for some unitary transformations C, P, and T satisfying = +T, = —C and 
= +P, acting on a Hermitian matrix Ai. For CI representations of the form 

(2.12), we identify C = a2, P = cr^, and T = —iai, where 0"^ denote the Pauh matrices 

that act on the N x N block matrices of Ai. 

The eigenvalue spectrum of Ai is reminiscent of the Wigner semicircle law, but the 

2N eigenvalues of A4 come in opposite-sign pairs ±1/^, with < z^i < . . . < z^at. Taking 

Q{0, a) to be a normal distribution, the joint probability density of the positive 

eigenvalues is given by 

a=l a<b a=l / 

In the Coulomb gas picture, the term "^^=1 ^a. encodes a repulsive force between each 
mirror pair of eigenvalues, ±z/a. In the eigenvalue spectrum this leads to a cleft at the 
origin, which is a subleading 1/A^ effect. The eigenvalue spectrum is shown in Figure 
2. 

The Altland-Zirnbauer CI ensemble has played a prominent role in the study of 
critical points in supergravity [19]: the critical point equation daVp = may be written 
as an eigenvalue equation for a matrix M. of the form (2.12),^ where Zab = Zi,a = 
T>aDbW is the fermion mass matrix. From studies of type IIB flux vacua [20], it is 
well-motivated to take the distribution (2.13) of the entries of Zab-, i-e. ^2(0, a), to be 

•^Our equation (2.12) differs from the definition oi M. in [19] by the phase of W , cf. equation (4.10), 
but the distinction is unimportant here. 
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uniform with support up to the flux scale. However, by universality only the first few 
moments of the distribution Q{0,a) are important at large N [13-17], and the exact 
distribution of the entries of Zab is therefore of little significance. 

With this identification of Zab, we immediately see that the positive eigenvalues Ua 
of M are identical to the diagonal entries of the Takagi factorization of Z presented 
in §2.1. Thus, the scalar mass spectrum in a supersymmetric vacuum can be expressed 
in terms of the eigenvalue spectrum of the Altland-Zirnbauer CI ensemble. We will 
find it useful to relate the latter to the spectrum of a certain Wishart ensemble, as we 
now explain. 

2.2.3 The relation between the CI ensemble and the Wishart ensemble 

To conclude this review, we will point out a formal equivalence between the squares 
of the eigenvalues of the CI distribution, and the eigenvalues of an almost square real 
Wishart matrix. The joint probability measure of the CI ensemble (2.17) is 

N N N 

f{v,, dv^ = 2-^C exp(-— Y.'^a + Y. l^l'^' - n ■ (2-18) 

a=l a<b a=l 

Defining = z/^ and C = 2~^C, this takes the form 

^ N N N 

\)l[dfia, (2.19) 

a=l a<b a=l 

which is the joint probability measure for the Wishart ensemble (cf. equation (2.9)) 
with P = 1 and ^ = 0. Since = M - N + 1 - 2/(5, this implies that M = + 1, so 
that the distribution of the squares of the eigenvalues v of the CI ensemble is equivalent 
to the distribution of the eigenvalues /i of the ensemble of real, almost square Wishart 
matrices. 

A further clarification is appropriate. The matrix ZZ is the Hermitian square of the 
square matrix Z, so one might incorrectly suppose that with the entries of Z distributed 
according to (2.13), the ensemble of matrices ZZ corresponds to the complex Wishart 
ensemble with M = N. However, Z is symmetric, while the matrix A appearing in 
the defining relation (2.8) has no special symmetry properties. What we have just seen 
is that in fact the ensemble of matrices ZZ corresponds to the real Wishart ensemble 
with M = A^ + 1. 

3 The Spectrum of the Hessian Matrix 

In the preceding section we showed that the scalar mass spectrum in a supersymmetric 
vacuum is completely determined (for fixed vev \W\) by the eigenvalue spectrum of 
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the Altland-Zirnbauer CI ensemble, which in turn can be obtained from the spectrum 
of the real Wishart ensemble with M = N + 1. In this section we will use existing 
results for the Wishart ensemble to obtain an analytic expression for the scalar mass 
spectrum. 

In §3.1 we will show that an AdS vacuum with ^ 1 and W ~ Zab has a large 
number of BF-allowed tachyons, which may be problematic for model building. In §3.2 
we compute the fraction of vacua in which there are no tachyons, and we obtain the 
mass spectrum in these rare vacua. 

3.1 The spectrum of Hz from the Marcenko-Pastur law 

From (2.6) we note that the spectrum of Tiz interpolates between that of Marcenko- 
Pastur (for small vevs of \W\), and that of the CI ensemble (for \W\ — ?■ oo). More 
precisely, this interpolation is controlled by the dimensionless ratio where 



susy 



is the supersymmetric mass scale, i.e. the scale of the entries of Zab, and^ ^7^3/2 



\W\/Mpi. To be fully explicit, we will temporarily reinstate factors of Mpi, and we 
factor out a scale m^usj/ from the tensor Zab as 

Zab = TTT'Susy Zab, (3-1) 

where rrisusy is determined so that^ Zab = £ ^(0, I/a/ZV), for a ^ b, and Zaa G 
f2(0, ^J2/N). Finally, we introduce the notation 

\W\ = m3/2 = rrisusy W Mli , (3.2) 

for a dimensionless parameter W. The eigenvalues of the matrix M = Ai /rrisusy are 
correspondingly denoted by A, and we note that the spectrum of each 2x2 block of 
Tiz can be written as (by here and henceforth suppressing indices) 

UJ± = ml^sy = ^lusy (^^ ± l^^l^) , (3-3) 

with A > 0. We will now obtain the distribution of 0;^ for arbitrary values of \yV\- 
The probability that Cj± is smaller than a given value a is given by 



P (cl;± < a G Pi) = P ± |I?|A < a) = P 1^ 



A±M 



P(A<T^ + v/a + T A<T^-v/a+T) , (3.4) 



^Recall that in each vacuum the expectation value of the Kahler potential has been set to zero. 
Choosing the standard deviation of Zah to be proportional to \/^fN yields mass spectra that are 
A^-independent to leading order in 1/iV; the standard deviation of the Z^b can then be changed by 
rescaling rususy 
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where the domains of support for the probabihty density functions of the two branches 
of Hz are given by T'+ = [0, oo) and = [— |PF p/4, cxd), respectively. Due to the 
positivity of A, the last term of (3.4) does not contribute to the probability of the 
positive definite branch, and it only contributes to the a)_ branch for a < 0. 

Equation (3.4) demonstrates how the probability density of the eigenvalues of the 
supersymmetric Hessian matrix is completely specified by the distribution of A, which 
is given by the Altland-Zirnbauer CI ensemble. But as shown in §2.2.3, the joint 
probability density of the squares of the eigenvalues in the CI ensemble is equivalent 
to that of the eigenvalues of 'almost square' (M = N + 1), real Wishart matrices, as 
in equation (2.19). Thus, the probability distribution of the eigenvalues of the Hessian 
matrix may be written in terms of the Wishart eigenvalue fi = as 



where /3_ = ^^y- 



Picd+ <d eV+) = P (fi< /3_ 



^ -r a 

Pi 



-|I?|^a + ^, and 



P{uj- < a eV^) = P {jl < (3+^ - P (^fi < (3- a < , 

J — h a + \ W\\J~h^-^^-. The probability density of the t 
l-Lz is then given by 

_ dP{Cj+ <a) _ dP{fi </?-)_ ^ ^ , 
da da dP- da ' 



with /3+ = ^ + a + |1? 
igenvalues of 



ei; 



(3.5) 
(3.6) 

two branches of 



for 



_ dP{li < /3_) 
da 

deV^ 



Pu^+id) 
and 



da 

dP{Lo_<d) d(3+ dP{ji < (3+) d(3_ dP{fi < (3_\d < 0) 

Puj^[a) = j2 = -t:^ ^ j2 

da da dB< da 



(3.7) 



= ^ PMP0+) - 0(-a) ^ PMP0-) , yo.o) 

for a G T>^, where Pmp{P±) is given by (2.10) for M = + 1. The spectrum of the 
ensemble of Tiz matrices expressed in terms of a (i.e. in units of the supersymmetric 
"^""s scale) is then given by 

1 . 



da dp+ 



dp_ 



(3.J 



mass I 



PHzi^) = 



1 



1 / .A 
= 2 I ^ P'^'^^- 



da 



- 0(-a) ^ PMP0-] 



1 



\ + . I^L \ PMP0- 




tj;j,„p(;3_)J,(3.9) 
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which has support in the domain a G = [— |l^p/4, oo). Upon including the shift 
given in equation (2.5), the probabihty density of the scalar mass spectrum is 

PH{m') = pn,{m' + 2\W\^) , (3.10) 

where py_^ is given in equation (3.9). This is one of our principal results. 

Some features of this distribution of supersymmetric masses are worth extra atten- 
tion. First, we have shown that the spectrum of the Hessian matrix is given in terms 
of a certain Marcenko-Pastur law, and features of the Wishart spectrum give rise to 
characteristic features in the spectrum of l-iz and of the Hessian "H. 

Expanding the edge positions rj± given in equation (2.11) to leading order in 1/iV, 
Pmp{.p) of equation (2.10) for M = 1 is nonvanishing only between the 'hard edge' 
at /t = and the 'soft edge' at /i = 4. Since Cj± = (i± \ W\^/]l, this means that the 
positive branch of l-iz has support for 

Q<u+<Aml^^y + 2\W\m,usy. (3.11) 

to leading order in 1/A^. The domain of support of the negative branch depends on 
\W\: for > Arrisusy, the probability density function is nonvanishing for 

4mL^ - 2\W\msusy < oo- < , (3.12) 

while for \W\ < 4 rrisusy, the spectrum develops a hard edge at U- = — |iyp/4, and has 
support in the range 

- ^ < < Max(0,4mL, - 2|Ty|m,„,^) . (3.13) 

In particular, for > 2 rususy the domains of support of the two branches do not 
overlap. 

The spectrum of the Hessian is given by a simple translation by — 2|14^p of the 
spectrum of Tiz, and by equation (2.7), the hard lower edge of equation (3.13) cor- 
responds to the BF bound for eigenvalues of the Hessian. Note that after the shift 
the negative branch has only negative support if \W\ > rrisusy, and for \W\ > 2 rrisusy 
all masses are tachyonic. Also, by equation (3.9), for \W\ < Arrisusy the probability 
density pi-i{m?) has a square root divergence at the BF bound, which leads to a sub- 
stantial amount of tachyons with masses close to the bound. We have plotted the 
typical mass spectrum for a supersymmetric AdS vacuum for the six different values 

= '^^rrisusy) \m susy ^rn susy, '^rrisusy) 5 rrisusy: 20 TH-s^jsj, in Figure 3. One can see that 
the shape of the spectrum interpolates between the Wishart spectrum (cf. Figure 1) 
for small \W\ and the Altland-Zirnbauer CI spectrum (cf. Figure 2) for large \W\. 
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Figure 3. The mass spectra, in units of rri^^^y, of generic supersymmetric AdS vacua with 
= 100 complex fields. The purple (darker) regions correspond to the contribution from the 
positive branch of (3.3), and the blue (lighter) regions correspond to the contribution from the 
negative branch. The black vertical line on the left side is at the BF bound w-gp = — 
(not shown for plot 3(f)). 
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Using equations (3.5) and (3.6), we find that the total fraction of scalar fields that 
are tachyonic can be expressed in terms of the Wishart eigenvalue fi as 



ftachyons = ^ (A < M^f) + ^ (A < 1^^!') ] , (3.14) 

which e.g. for \W\ = nisusy is around 80%. 

The cleft in the Altland-Zirnbauer CJ spectrum translates into a similar cleft in the 
supersymmetric mass spectrum at m? = —2\W\'^. As in the case of the CJ ensemble, 
the cleft is a sub leading effect in 1/A^, with a width that is inversely proportional to 
N. Moreover, as the probability distribution of goes to zero at the right edge of 
the Wishart spectrum at a = C+ = 4 + 0{1/N), there is a kink in the supersymmetric 
mass spectrum at = '^m'l^^y — 2\W\{mausy + \ We have verified these features 
through explicit numerical simulations, as detailed in the appendix. 

3.2 Fluctuations to positivity 

As we have just shown, typical supersymmetric vacua in supergravity have BF-allowed 
tachyons. In anti de Sitter space, these tachyons do not destabilize the vacuum, but in 
order to obtain a background with a positive cosmological constant the AdS solution 
needs to be uplifted. While it is not impossible that a suitably constructed uplift 
potential may stabilize several of the tachyons, requiring that the resulting de Sitter 
critical point is tachyon free is a significant restriction on string theory model building. 
It is therefore interesting to consider the supersymmetric AdS vacua that are free of 
tachyons, or more generally have masses only above a certain bound. This subset 
of solutions would appear more likely to result in metastable de Sitter vacua after 
uplifting. 

In this section we compute the fraction of supersymmetric vacua with masses above 
any given positive bound by mapping the problem of the positivity of the eigenvalues 
of the Hessian to that of fluctuations of the corresponding real, almost square Wishart 
matrix [21, 22]. We present an analytic formula for the probability of such fluctuations 
as a function of the minimal mass squared, \W\, rrisusy, and N, and we derive the 
corresponding mass spectrum. 

3.2.1 The probability of positivity 

As we have discussed above, the scalar mass spectrum of the ensemble of supersym- 
metric vacua is determined in terms of the eigenvalues of a real Wishart matrix with 
M = + 1. At large A^, the eigenvalue spectrum of the Wishart matrix is given by the 
Marcenko-Pastur law, which is bounded by /t G [r7_,?7+], as in equation (2.10). How- 
ever, while the Marcenko-Pastur law determines the typical spectrum of eigenvalues. 
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more rare configurations exist. Here, we are interested in fluctuations of tlie smallest 
Wishart eigenvalue fimin such that /imin > 4|iyp. In this case the smallest eigenvalue 
of the Hessian T-L is given by 

"^min(Amin) = "^Lsy " (^min " \W\^/ flrain - 2\W\^^ > 0, (3.15) 

and the mass spectrum is tachyon free. Note that for \W\ > rrisusy these fluctuations 

require that all Wishart eigenvalues fluctuate past the generic right edge of the spectrum 

at r]+=A + 0{l/N). 

The probability of a fluctuation to positivity is therefore given by the probability 

that the smallest eigenvalue of the corresponding real Wishart matrix, with M = 
+ 1, fluctuates to /tmin > 4|W^p. In the case of real Wishart matrices with M = 
+ 1 for which the entries of A appearing in (2.8) are normally distributed, Edelman 

has obtained the exact distribution of the smallest eigenvalue, making no large 

approximation [7]: 

PA.J^) = N'xliN'x) = ^e-'2^'\ (3.16) 

Universality very plausibly leads to compatible results for more general distributions 
of the entries of A, once is sufficiently large. 

The probability that the smallest eigenvalue fimm is larger than some value ( is 
then given by 

POO 

Pifimin >0= dxpf,_{x) = . (3.17) 



_ 1 

Thus, the probability that a random, supersymmetric AdS vacuum has no tachyons is 
P(m^ > 0) = exp (-2Ar2||^n = e^^{-2N'^\W\^ /ml^^y) . (3.18) 

This is one of our main results. Unless \W\ < 0{1/N), fluctuations to positivity are 
extremely unlikely at large A^. 

3.2.2 The fluctuated spectrum 

The eigenvalue spectrum that results from a fluctuation to positivity can be computed 
using the Coulomb gas method as discussed in e.g. [23, 24]. In [22], Katzav and Perez 
Castillo studied Wishart matrices for which all eigenvalues satisfy fi > ( > rj^, flnding 
the eigenvalue spectrum (for a = 1/ \/N and M = A^ + 1) 

Pmp(A) = ^ ; , (3.19) 
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with 



9(0 



4(2 + Ar(4 + C)) 
3A^ 



sm 



TT 1 

— I — arccot 
6 3 



27NC 



{2 + N{A + C))^ -27NC 



(3.20) 



At the left, hard edge ( of the fluctuated Wishart spectrum, the spectrum (3.19) 
has a characteristic square root divergence. The right edge of the fluctuated Wishart 
spectrum is at /tmax = g{(), which for ( = 4|iyp is at /imax = 4(1 + \W\'^) + 0{1/N). 

As shown in §3.1, the spectrum of Hz is given in terms of the (now fluctuated) 
spectrum of the corresponding Wishart matrix. By conditioning on fi > C = 4|iyp, 
the second term in (3.6) vanishes, so that the fluctuated spectrum of Tiz is given by 



\w\ 



A/4a+|T?p 



PmUM + 0(a) 1 



\w\ 



PmpW- 



(3.21) 



where again 



a±\W\\/a+^-^, andplTpip) is given in (3.19), (3.20). From 



2 i-"^|"|y--i 4 

equation (2.5), we note that upon including the shift, the probability density of the 
fluctuated scalar mass spectrum is 



Pn'ira') 



p^^l{m' + 2\W\ 



The two branches of eigenvalues of "H have support for 

C±\W\^- 2\W\^ < < 4 + 2\W\ (\W\ ± ^1 + \W\'^\ 

^susy \ / 



(3.22) 



(3.23) 



Note that for C = 4|iy|2, this means that the branches do not overlap when \ W\ > -y=. 
The correlation between the two branches leads to two distinct peaks in the mass 
spectrum, at each of which the spectrum exhibits a square-root divergence. In Figure 
4 we have plotted the fluctuated spectrum for ( = 4|W^p and for two different values 
of \W\. 

3.3 The distribution of tachyon-free vacua 

Let us conclude this section by considering the full distribution of supersymmetric 
vacua, and the constrained distribution of tachyon-free supersymmetric vacua, in ran- 
dom supergravity. The number of vacua can be computed as [20] 



N.ac. = J dV[W, F,Z,.. .] 5(2^) (F) I det BaFbI 



(3.24) 



where the capital indices A, B run over both holomorphic and antiholomorphic indices. 
The measure dV specifies the prior distribution of the fully covariant Taylor coefficients 
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(a) = I msusy (b) \W\= nisusy 

Figure 4. The mass spectra, in units of m'^^gy, of supersymmetric AdS vacua with = 100 
complex fields that have fluctuated to positivity. The purple (darker) regions correspond to 
the contribution from the positive branch of (3.3), and the blue (lighter) regions correspond to 
the contribution from the negative branch. Blue curve: analytic result from equation (3.22). 
Black line: equilibrium eigenvalue distribution from Metropolis simulation (see Appendix A). 



of the superpotential, and as in the rest of this paper, we will make the assumption that 
the integration variables are independent and identically distributed, while allowing for 
different standard deviations for \W\, \Fa\ and \Zab\- In particular, this means that the 
prior measure is separable and, with some abuse of notation, we write dV[W, F, Z] = 
dV[W]dV[F]dV[Z]. With these assumptions we have 



N^ac. =Cn j dV[W] j dV[F] j dV[Z] 5(2A^)(F) \ det OaF, 
= Cn j dV[W] J dV[F] J dV[Z] 5(2^^)(F) 



. (3.25) 



The integral over Zab can be written as an integral over the eigenvalues and eigenvectors 
of the matrix Ji4 in equation (2.12). By applying the 2A^ x 2A^ unitary transformation 
of equation (2.3), it is easy to see that the above integral simplifies to 

N.ac. = CN J dV[W] j dV[F] /"dP[Ai,...,A;v] ni^'-l^l'l ' (3-26) 

a=l 

where the integration over eigenvectors of M. is suppressed. The Jacobian factor 
l-^a ~ I^PI ensures that each vacuum counts with unit weight in the integral 
by canceling the corresponding inverse factor arising from the delta function. 
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The frequency of tachyon-free supersymmetric vacua is given by 



_ f dV[W] fdV[F] f dV[\i, . . . , Xn] S^^^\F) Ua=i\^l-\W\^\ e-^^'l^l'/"^-^ 
" / dV[W] JdV[F] /rfP[Ai,...,A^] 6(^^)iF) nf=i|A^-|W^n ' 

where nisusy is determined as in the discussion around equation (3.1). 

Of particular interest is the region in which Ai > or equivalently rrtsusy ^ 
A^|iy|. In this region there is no exponential suppression of tachyon-free vacua, and 
in fact tachyon-free vacua are quite common. For dV[W] being roughly uniform in the 
complex W plane, a fraction l/N"^ of all supersymmetric vacua, and most tachyon-free 
supersymmetric vacua, occur in this region. The spectrum in this more stable region 
consists of two overlapping branches approximately following the Marcenko-Pastur law, 
as illustrated in Figure 3(a). However, the absence of BF-allowed tachyons in a super- 
sjTumetric solution does not guarantee the existence of an uplift to a metastable de 
Sitter vacuum, as we will now discuss. 



4 Instabilities of Uplifted Vacua 

Our results so far concern the prevalence of tachyons in supersymmetric vacua, but 
the question of primary physical interest is the likelihood of tachyonic instabilities in a 
realistic vacuum with positive cosmological constant. The 'uplift' paradigm advanced 
in [1] suggests dividing the problem into two parts: ascertaining the mass spectrum in 
a supersymmetric vacuum, and incorporating stabilizing or destabilizing effects of the 
supersymmetry-breaking energy. In §3 we have given a comprehensive treatment of the 
first point for A/" = 1 supergravities that are 'random' in the sense of [2], and we now 
turn to quantifying the effects of uplifting. 

Instabilities of random critical points in various models of a supergravity landscape 
have been studied in a number of works, which we briefly review for completeness. 
The original analysis of [25] modeled the mass matrix at a critical point as the sum 
of a real Wigner matrix and a positive definite diagonal matrix, and quantified the 
probability of instability in terms of the size of the Wigner component. More recently, 
[8] worked with a similar model of the mass matrix, motivated by constructions in 
type IIA string theory. The stability of de Sitter vacua resulting from spontaneous 
breaking of supersymmetry by an F-term was the subject of [2, 19]. In those works the 
supersymmetry breaking was taken to result from the same W and K that gave rise to 
the moduli potential, rather than through an uplifting due to supersymmetry breaking 
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in another sector. In general, explicit uplifting from a supersymmetric solution to a 
de Sitter solution can arise by integrating out the supersymmetry breaking sector of a 
theory with spontaneously broken supersymmetry, but not all conceivable uplifts need 
to be of this form. 

Given the mass spectrum in a supersymmetric vacuum, the crucial question is 
whether the source of uplifting introduces new instabilities, cures existing instabilities, 
or leaves the mass matrix unchanged. A very simple and presumably unrealistic picture 
that is often invoked in the literature is a 'rigid' uplift, in which the cosmological 
constant is increased but the mass matrix is unmodified. In this case, our counting 
of tachyons in supersymmetric vacua translates directly to the de Sitter critical points 
resulting from uplifting. In particular, in view of equation (3.18), the probability that 
a rigidly-uplifted de Sitter configuration is metastable is 

P{m^ > 0) = exp (-2iV2|iy|VmL,) , (4-1) 
so that metastability is extremely rare for \W\ ^ nisusy/N. 
4.1 Wigner uplift 

While modeling the uplift as rigid serves as a simple starting point, models in which 
the uplift potential has a nontrivial field dependence are more general and appear more 
likely to arise from supersymmetry breaking in string compactifications. For a general 
function Vup of N complex scalar fields, the contribution to the Hessian is given by 

[dl^updlv^pj ' ^"-'^ 

where d'^jVup = (<9|feKip)* is an x Hermitian matrix and S^^Vup is an A^ x A^ 
complex symmetric matrix. We will take the entries of d'^jVup to be distributed as the 
i.i.d. entries of an A^ x A^ Hermitian matrix with the natural scale V^p ~ 3|iyp, and the 
entries of dlf^Vup to be distributed as the i.i.d. entries of an A^ x A^ complex symmetric 
matrix of the same scale. 

We note that Hup is 'time-reversal' symmetric, 

TnupT^ = n:p, (4.3) 

under the unitary time-reversal operator T = +T~^ = ai acting on the N x N block 
matrices. Since this is the only symmetry of Hup, we identify Tiup as a representation 
of the A J symmetry class [6, 26], more commonly known as the Gaussian orthogonal 
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ensemble. Thus, although we have taken Vup to be a function of complex scalars, Hup 
is in the symmetry class^ of a real symmetric matrix, just as in the model of [25]. 

Supersymmetric vacua in which irisusy ^ A^|W^| have many tachyons in AdS, and 
as a Wigner uplift potential of the form (4.2) serves to disperse the supersymmetric 
spectrum, these vacua will predominantly uplift to unstable critical points in Minkowski 
or de Sitter space. The probability of obtaining special configurations that enjoy a 
larger likelihood of uphfting to a metastable vacuum was discussed in §3.2, and here we 
simply note that in order to produce a vacuum once the uplift potential is taken into 
account, the spectrum of the AdS vacuum typically has to fluctuate far enough so that 
the left edge of the spectrum is well above zero, which is extremely unlikely in general. 

A much more interesting region occurs for \W\ < in which case the AdS 

vacuum is typically tachyon free, and the uplift potential will only produce a minor 
perturbation in the spectrum. We find it convenient to parametrize the standard de- 
viation of the entries of 1-iup as ^Uup = C2m^„g^/A^^, for a dimensionless constant C2, 
and we parametrize \W\ as \W\ = Ci/N. Note that the natural scale Vup ~ 3|W^p 
corresponds to Ci ~ C2. Specializing to a basis in which the supersymmetric Hessian 
is diagonal, "H = diag(m^), it is easy to compute the eigenvalues of the total Hessian, 
'Htot = T-L + Tiup, to second order in perturbation theory: 

= ml + {Hup)u - E tII?^ + • • • , (4.4) 

where the ellipsis includes terms of higher order in l-Lup- The typical size of the smallest 
eigenvalue can be estimated from (3.15), using the fact that from (3.16), the mean of 
the smallest eigenvalue of the associated Wishart matrix is 2/N'^: 

2 

777 

^^{2 -a), (4.5) 

where a = v^ci + 2c^ + a/2c2, and the quadratic corrections have been neglected. Thus, 
to leading order a typical de Sitter critical point arising from uplifting remains stable 
if ci, C2 < 1. 

It remains to check that higher order corrections are sub leading. For small \ W\, the 
eigenvalues mf are nearly doubly degenerate, and the two smallest eigenvalues are split 

^In our numerical treatment we have simulated matrices of the form (4.2) and found perfect agree- 
ment with the Gaussian orthogonal ensemble, strongly suggesting that these matrices are in fact in 
the same universality class. 
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Figure 5. The probability of the absence of tachyons in H-tot = T~L + 'Hup for \W\ = ci/N, 
versus ci for A^=10 and C2 = 2\/3ci, \/3ci, ci, (from left to right). 

by {ml - ml) = 2^/2\W\ml^,.y/N = 2V2cim1^^y/N^. For ci < 1 this is the smallest 
difference of eigenvalues in the spectrum, and the expectation value of the quadratic 
correction to the smallest eigenvalue can be approximated by (cf. (4.4)) 



so that higher order corrections are indeed subleading for ci ~ C2. 

We have verified these analytic result through extensive numerical simulations. 
The probability of the absence of tachyons in the mass spectrum of T-Ltot = 'H + l-Lup foi' 
different values of Ci and C2 is shown in Figure 5. 

4.2 D-terms 

Although it is plausible that many sources of supersymmetry breaking and uplifting 
will introduce instabilities, as described in the preceding section, we now point out that 
the presence of D-terms can contribute a stabilizing effect. A comprehensive treatment 
of D-terms in random supergravity is beyond the scope of this work, and we will content 
ourselves with a few observations. 

For an A/" = 1 supergravity with the product gauge group G = Gi x G2 x • • • x Cn, 
the scalar potential contains the additional term 



\ [rnf -ml) / ^ 2^ci/m ^ Ci ' 




(4.6) 



Vd 



(4.7) 
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where fij denotes the holomorphic gauge kinetic function, and Di is the D-term for the 
gauge group Gi, which (for field configurations with W 0) is given by [27] 



i 1 



a 



(4.8) 



where e*X" is the variation of the field 0" under an infinitesimal gauge transformation 
A* — )■ A* + de\ Here, denotes the components of the Killing vector of the corre- 
sponding isometry of the Kahler geometry, which for linearly realized gauge symmetries 
is given by X" = —iTi(j)°', where we have suppressed the gauge indices for the generator 
Ti and for 0". 

A supersymmetric solution satisfies DaW = 0, so that it follows from (4.8) that Vd 
vanishes at any supersymmetric critical point. For this reason D-terms cannot be used 
to uplift an F-fiat minimum. However, the D-term potential can change the masses 
of the scalar fields. In particular, since Vd is positive semidefinite and vanishes at 
a supersymmetric vacuum, it can only contribute positively to the eigenvalues of the 
Hessian, and may therefore decrease the likelihood of tachyons. However, the presence 
of D-terms restricts the contribution to the Hessian matrix arising from Vp as follows: 
since the D-terms are real, (daDi)* = daDi, which for supersymmetric solutions implies 
that^ ZabX^ = —Sg^i^WX^, where we have specialized to K^i = at the point in 
question. This consistency condition can be written as an eigenvalue equation for the 
2N component vectors Xf = (Xf , ±Xf )^ as 



has the same spectrum as the matrix A4 of equation (2.12), and was identified in 
[19] as the matrix governing the critical point equation of spontaneously broken F- 
term supergravity. Here 'dw denotes the phase of the superpotential, W = e*''^|Vr|. 
It follows from equation (2.6) that the contribution to the Hessian from the F-term 
potential will have eigenvalues and — 2|iyp in the directions Xf. 

Using DaW = and redefining the fields A* to set Refij = 6ij at the point in 
question, we find the D-term contribution to the Hessian in a supersymmetric vacuum: 



MXf = t\W\X^, 



(4.9) 



where 




(4.10) 



Xia X'^^ —Xia X^^ 
—Xia Xia X\ 




^This relation may also be obtained directly from the Killing equation. 
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where the indices a, b are lowered by the metric and the index i is raised by the gauge 
kinetic function. The D-term contribution to the Hessian for any given Gi is thus of 
rank one, with a nonvanishing contribution only in the directions specified by . 

We conclude that for n N, D-terms do not substantially alter the results de- 
scribed in §4.1, but do lead to the existence of special directions along which the 
eigenvalues of the Hessian are given by and — 2|iyp + 

In cases where n > N, contributions from the D-terms may have a substantial 
impact on the mass spectrum and the number of tachyons. In order for the D-terms 
to affect the stability of the Hessian arising from the F-term potential, a substantial 
fraction of the scalar fields have to carry charges under the gauge groups Gi. Under- 
standing the extent to which this condition is met in well-motivated compactifications 
is an interesting problem that we will not address in this work. 

5 Conclusions 

We have computed the scalar mass spectrum in a supersymmetric vacuum of a general 
four-dimensional Af = 1 supergravity theory, with the Kahler potential and superpo- 
tential taken to be random functions of ^ 1 complex scalar fields. By relating 
the spectrum of the Hessian matrix to the spectrum of an associated Wishart matrix, 
we showed that a fraction / = exp{—2N^\W\'^ /ml^^y) of AdS vacua are tachyon-free. 
Then, using results from the Coulomb gas formulation of random matrix theory, we 
obtained the scalar mass spectrum resulting from a fluctuation to positivity. We per- 
formed extensive numerical cross-checks of our analytic expressions for the fluctuation 
probability and mass spectrum, with perfect agreement. 

A clear implication of our results is that uplifting a supersymmetric vacuum with 
\W\ > rrisusy/N through the addition of positive-energy sources is overwhelmingly 
likely to lead to an unstable de Sitter critical point, not a metastable de Sitter vacuum. 
The likelihood of uplifting to a de Sitter vacuum can be increased by making the 
supersymmetric mass scale rrisusy large compared to \ W\. For rrisusy N\W\ and A^ ^ 
1, the probability that uplifting a supersymmetric AdS vacuum leads to a metastable 
de Sitter vacuum can be vastly larger than the probability that a general de Sitter 
critical point with spontaneously broken supersymmetry is metastable, cf. [2]. 

Our results are applicable when the superpotential and Kahler potential are accu- 
rately described as random functions of A^ complex fields. Determining the prevalence 
of tachyons in AdS vacua in more general theories with multiple sectors or special 
structures in K and W is an interesting problem for the future. 
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A Numerical results 

In §3 we obtained an analytic expression for the supersymmetric mass spectrum l-iz 
and computed the probability of a fluctuation to positivity. The simplest and most 
direct check of these results, the direct diagonalization of a large sample of Wishart 
matrices, is a computationally intensive task: most of the parameter space of interest 
corresponds to large fluctuations that are extremely rare (i.e. |iy| ~ 1 and A ^ 1 in 
(3.18)). A more efficient approach utilizes Monte Carlo methods to obtain both the 
fluctuation probability and the spectrum of fluctuated eigenvalues. In the following, 
we present two numerical techniques to obtain these quantities. 

To begin, we will review a few elements of the Coulomb gas formulation of random 
matrix theory [26, 28], which underlies the analytic results presented in §3.2, but also 
serves as an efficient starting point for simulations [24]. 

Recall that the joint probability density function of the eigenvalues of a Wishart 
matrix is given by (2.9). The integral over this joint pdf can be interpreted as the 
partition function of a gas of interacting particles undergoing Brownian motion in one- 
dimensional space at finite temperature. The eigenvalues are then interpreted as the 
positions of particles that are confined by a linear background potential and repel each 
other with a force that is inversely proportional to the distance between any pair of 
particles.^ The Hamiltonian of this N particle system is then given by 

^ N N N 

H = -^^lla -'2^^M^^a- l^b\- i^^'^^J^a, (A.l) 

a=l a<b a=l 

where we impose the constraint > 0, which corresponds to a hard wall at /i = 0. 
The temperature of the gas is given by 2//3 — recall that /3 = 1, 2 for real and complex 
Wishart matrices, respectively. 

^This force may be interpreted as a two-dimensional Coulomb interaction between equally charged 
particles moving in one-dimensional space, hence the name of the formalism. 
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The Coulomb gas formalism suggests that techniques and intuition from statistical 
mechanics can be useful in the study of random matrices. In fact, by rewriting the 
partition function in terms of the empirical eigenvalue density p(/i) = ^ J2iLi f^i)^ 
and taking the continuum limit, the partition function can be written as a path integral 
over the scalar field p, and saddle-point evaluation can be used to find the equilibrium 
distribution (2.10). 

Moreover, saddle-point evaluation can be used to find the probability density of the 
eigenvalues of constrained configurations, such as those in which the smallest eigenvalue 
is larger than some cutoff (, as well as the probability of obtaining such a configuration. 
In §3.2, we used solutions of the fluctuated Wishart spectrum obtained in [22] to find the 
spectrum of fluctuated eigenvalues of the Hessian matrix analytically. In this appendix, 
we confirm these findings by numerical simulation of the A^-particle Coulomb gas. 

The probability of obtaining a fluctuated solution is given by 

^(/^min > C) = , (A.2) 

Zn{-oo) 

where the constrained partition function ^Ar(C) is given by 

/oo POO POO / 1^ \ 

rfpi j dn2... J duN exp y--H{n)j . (A.3) 

As in [24], the integral in (A.3) can be computed numerically via Monte Carlo sampling 
by considering the average (■)^, where the Pa > C ^i^^e randomly chosen from a uniform 
distribution. This leads to the probability of the smallest eigenvalue 

where the prefactor captures the average spacing between the sampling points. We used 
6 X 10^'' Monte Carlo samplings to evaluate the fluctuation probability numerically and 
obtained excellent agreement with the analytical result, as shown in Figure 6. 

While the method of evaluating the partition function directly via Monte Carlo 
integration captures the fluctuation probability, it does not provide the probability 
density function of the eigenvalues. In particular, we are interested in the pdf of a 
fluctuated spectrum, as computed in (3.19). However, these fluctuations are extremely 
rare for the parameters of interest, and are essentially inaccessible by direct diagonal- 
ization of matrices. Instead, we will take advantage of the Coulomb gas picture and 
numerically simulate the equilibrium distribution. 

A popular algorithm to thermalize classical particles in an arbitrary potential is 
the Metropolis algorithm [29]. This algorithm shifts one particle at a time to obtain 
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Figure 6. The probability of a fluctuation to positivity of Hz for \ W\ = 1, versus the number 
of complex scalar fields N . The curve is the analytic result given in equation (3.18), and the 
dots correspond to 6 x 10^*^ Monte Carlo samplings with 3<t error bars. 

the equilibrium distribution in the following way: (1) select one particle, (2) obtain the 
energy cost for performing a random step, (3) if the energy cost is negative, keep the 
step. Otherwise, keep the step with probability e~^^-^/^. In order to avoid artifacts 
from the simulation, the size of each step should be small compared to the typical 
particle spacing. When starting with an arbitrary particle distribution there is a typ- 
ical number of steps ^ithermaiize required after which the initial correlations decay and 
the gas reaches thermal equilibrium. In our simulations we chose a fixed number of 
steps nstcps ^ ^thcrmaiizc to avoid correlations from the initial configuration but still 
observe any thermal fluctuations. While implementing the Metropolis algorithm for 
the Hamiltonian in (A.l) leads to a generic random distribution of eigenvalues, we can 
further constrain the Hamiltonian by introducing a hard wall at = allowing us 
to access the eigenvalue distribution for large fluctuations. The result from Metropolis 
simulation, along with the analytic result, appears in Figure 4. 
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